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ABSTRACT 

This  paper  investigates  the  computational  analysis  of  unsteady  flow  of  blood  and  heat  transfer  through  a  stenosed 
artery  in  a  third  grade  fluid  model  with  slip  conditions.  Incorporated  into  the  model  is  an  externally  applied  magnetic 
field.  The  solutions  of  the  unsteady  non-linear  dimensionless  momentum  and  energy  equations  are  obtained  using 
Galerkin’s  weighted  residual  and  Fourth  order  Runge-Kutta  methods.  Effects  of  slip  velocity,  magnetic  field,  shear 
thinning,  shear  thickening  and  other  parameters  on  the  flow  and  heat  transfer  characteristics  are  presented  graphically 
and  discussed. 
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INTRODUCITON 

The  study  of  blood  flow  in  the  cardiovascular  system  has  generated  a  lot  of  interest,  hence,  the  huge  amount  of  literatures 
are  studied  on  the  subjects.  Sanjeev  and  Chandrashkhar  [1]  considered  mathematical  model  of  power  law  fluid  with  an 
application  of  blood  flow  through  an  artery  with  stenosis.  Their  results  showed  that,  increases  in  the  size  of  the  stenosis 
lead  to  increase  in  the  pressure  drop  and  flux  but  decrease  the  resistance  to  flow  and  wall  shear  stress.  They  concluded  that 
the  non-Newtonian  behavior  of  blood  is  helpful  in  the  functioning  of  diseased  arterial  circulation.  The  results  concerning 
the  velocity  and  temperature  field  of  Newtonian  blood  flow  in  the  presence  of  applied  magnetic  field  was  presented  by 
Kumar  and  Diwakar  [2].  They  found  out  that  the  rate  of  heat  transfer  indicates  the  presence  of  magnetic  field  appreciably 
influenced  by  the  application  of  the  magnetic  field.  Dharmendra  [3]  presented  a  theoretical  investigation  to  examine  some 
of  the  significant  characteristics  of  three  layered  oscillatory  blood  flow  through  stenosed  arteries.  Their  results  revealed 
that  the  size  of  intermediate  and  peripheral  layers  reduces  in  the  expanded  region  with  the  increasing  viscosity  of  fluid  in 
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peripheral  layer,  whereas,  the  opposite  effect  is  observed  for  the  viscosity  of  fluid  in  the  intermediate  layer.  The  pulsatile 
flow  of  the  blood  as  two  fluid  models  with  the  suspension  of  the  erythrocyte  in  the  core  peripheral  of  plasma  as  a 
Newtonian  fluid  was  investigated  by  Devajyoti  and  Uday  [4].  From  their  finding,  increased  slip  velocity  lead  to  increase  in 
the  flow  velocity  and  flow  rate  but  decreased  the  effective  viscosity.  Effect  of  slip  velocity  on  blood  flow  through  a 
catheterized  artery  was  investigated  by  Narendra  et  al  [5].  They  assumed  that,  the  flowing  blood  is  represented  by  a 
Newtonian  fluid.  From  their  analysis,  they  concluded  that  the  slip  velocity  reduces  the  effective  viscosity  and  wall  shear 
stress. 


Furthermore,  Singh  et  al.  [6]  theoretically  investigated  the  influence  of  externally  imposed  periodic  body 
acceleration  of  the  flow  of  blood  through  a  time  dependent  stenosed  arterial  segment  by  taking  into  account  of  the  slip 
velocity  at  the  wall  of  the  artery.  The  effect  of  slip  velocity  on  blood  flow  through  an  arterial  tube  in  the  presence  of 
multiple  stenosis  was  studied  by  Arun  [7].  He  considered  the  effects  of  length  of  stenosis  and  shape  parameter  on 
resistance  to  flow  and  shear  stress.  He  observed  from  the  graphs  that,  the  parameters  have  small  variation  for  different 
values  of  stenosis  shape  parameter. 

Several  other  researchers  Bhatnagar  and  Strivastava  [8],  Amit  and  Shrivastaver  al  [11],  Raja  and  Varshney  [12]  to 
mention  but  a  few  have  in  one  way  or  the  other  modeled  and  studied  the  flow  of  blood  through  a  constricted  artery  under 
the  influence  of  magnetic  field  and  slip  conditions.  None  of  these  researchers  considered  the  unsteady  blood  flow  and  heat 
transfer  models  in  their  studies. 


However,  studies  on  unsteady  blood  flow  and  unsteady  heat  transfer  have  received  very  scanty  attention.  In 
unsteady  fluid  flow  or  heat  transfer,  the  conditions  at  any  point  of  the  flow  or  heat  transfer  are  time  dependent.  Thus,  the 
governing  equations  describing  the  fluid  flow  or  heat  transfer  become  cumbersome  to  handle  as  the  equations  now  have 
time  derivative.  Some  of  the  researchers  investigated  unsteady  blood  flow  model  without  considering  the  heat  flow 
includes  Geeta  and  Siddique  [13],  Mandal  [14],  Islam  [15],  Srikanth  et  al  [16]  and  Ikba let  al  [17]. 

Therefore,  in  this  study,  we  analyzed  the  unsteady  blood  flow  and  heat  transfer  through  a  stenosed  artery  in  a  third 
fluid  model  with  slip  velocity  at  the  constricted  artery.  The  presence  of  externally  applied  magnetic  field  is  also  taking  into 
consideration  in  this  study. 

Mathematical  Models 


The  momentum  equation  describing  the  unsteady  fluid  flow  and  unsteady  heat  transfer  as  obtained  by  Mohammed  [18]  are 
respectively  given  as 
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From  the  right  hand  side  of  equation  (2.1),  the  first  term  can  be  written  as 

p  fd2w  1  dw\  p  d  /  duA 

p\dr2  r  dr  )  prdr\  dr) 

while  that  of  equation  (2.2)  from  the  right  hand  side,  the  third  term  can  be  written  as 


(2.1) 

(2.2) 


(2.3) 
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K  fd2T  +  1  dr\  _  K  3  /  dr\ 
pCp  \dr2  r  dr)  rpCp  dr  V  dr) 


When  we  substituted  (2.3)  into  (2.1)  and  (2.4)  into  (2.2),  we  respectively  obtained 

dw  _  p  d  /  3w\  6 3  /3w\2  d2w  ^  2/33  /3w\2  +  a1  d2w  ^  a1  d3w  dP  o/3 q  ^ 

dt  prdr\  dr)  p  V dr)  dr2  pr  \dr)  pr  drdt  p  dr2dt  pdz  p 
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Since  we  employed  slip  velocity  in  the  constricted  artery  as  shown  in  Figure  1 ,  the  corresponding  slip  conditions 
to  equations  (2.5)  and  (2.6)  are  respectively  given  as 


w  =  ws  at 

r  =  R(z)  1 

dw  n 

—  =  0  at 

dr 

r  =  0  _ J 

(2.7) 

and 

T  =  T 

w 

at  r  =  R(z) 

?=" 

at  r  =  0 

(2.8) 

To  non-dimensionalize  equations  (2.5),  (2.6),  (2.7)  and  (2.8),  we  introduce  the  following  parameters  and 


variables. 


_  w 
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When  equation  (2.9)  is  substituted  into  (2.5)  and  (2.6),  after  simplifying  we  obtained  respectively 
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with  corresponding  dimensionless  slip  conditions  simplified  as 
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And 


dG  fdw\2  ,  fdw\4  ,  f  d2w\  (dw\  .  1  d  (  dG\ 
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as  the  dimensionless  energy  equation. 
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with  the  corresponding  dimensionless  slip  conditions  simplified  as 
9  =  0  at 


99  n  ♦ 
—  =  0  at 

dy 


y  =  Rb 
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(2.15) 


and  has  been  described  by  Young  [19]  and  Biswas  [20] 

«CE)  =  i__qi+£2f2E£l  fOT  |Z|<L 

Rq  2R0  l  L  V  1  '  — 

R0,  for  |z|>L  (2.16) 

Methods  of  Solution 

We  have  used  Galerkin  weighted  residual  method  for  getting  the  required  solutions  of  the  problems.  By  this  method,  to 
obtain  the  velocity  profile  of  the  unsteady  fluid  flow,  we  assume  a  solution  of  the  form 

w(y,  t)  =  a0(t)  +  ax(t)y  +  a2(t)y2 

Subjecting  (3.1)  to  the  slip  conditions  (2.12)  and  after  simplifying  yields 

Now, 
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Putting  (3.4)  into  (3.2)  to  obtain 
w(f,  t)  =  V01  r2  +  a0(t)(  1  -  r2)  +  a2(t)f2(  1  -  f2) 
For  convenience  sake,  we  drop  the  bar  and  write  (3.4)  as 
w(r,  t)  =  V01  r2  +  a0(t)(l  -  r2)  +  a2(t)r2(l  -  r2) 
The  residual  for  equation  (2.10)  can  be  written  as 


(3.4) 


(3.5) 


RRi(r,  a0(t),a2(t)) 


dw 

dt 


\  -  — - ■  Mr-  fO  -  n  (6  (fq 
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2  ^2^  I  2 

dr2  r 


/dw\3\  /ld2w  fl3w\ 

Car)  +  + 


(3.6) 


substituting  (3.5)  into  (3.6)  and  simplifying  we  obtain 

RRi(r,  a0(t),  a2(t))  =  -G  ~  ^](.4Voi  +  4a2(t)  -  4a0(t)  -  16 a2(t)r2)  +  d0(t)(l  -  r2)  +  d2(t)r2(l  -  r2)  + 
4121(a0(t)  -  a2(t)  +  4r2ri2(t))  -  (48 V0\r2  -  144V021a0(t)r2  +  144V^a2(t)r2  -  480V021a2(t)r4  +  144K01a^(t)r2  - 
40V01a0ta2tr2+960V01a0ta2tr4+144V01a22tr2-144V01a22tr44V01a22tr6+144a02a2tr2-480a02ta2tr4- 
144a0ta22tr2+960a0ta22tr4-  1344a0ta22tr6-  480a23tr4+1344a23tr6-  48a03tr2+48a23tr2-  1152a23tr8- 1 
28a23tr8+l  92V01a22tr6- 1 92a0ta22tr6+l  92a23tr6-  96V012a2tr4+l  92V01a0ta2tr4- 1 92V01a22tr4-  96a02 
ta2tr4+l  92a0ta22tr4-  96a23tr4+16V013r2-  48V012a0tr2+48V01a02tr2-  96V01a0ta2tr2+48V01a22tr2- 16 
a03t+48a02ta2tr2-  48a0ta22tr2+16a23tr2-Ml  V01r2+a0t-  a0tr2+a2tr2-  a2tr4 
(3.7) 


By  differentiating  (3.5)  with  respect  to  a0(t)  and  a2(t)  we  obtain  the  weight  functions  respectively  as 
Ui(r)  =  (1  —  r2)  (3.8) 

and 

U2(r )  =  r2(l  -  r2)  (3.9) 

We  obtained  the  following  systems  by  taking  the  orthogonality  of  the  residue  RR^dQ^t),  d2(t),  r)  with  respect  to 
the  weight  functions  t/1(r)and  U2(r ) 

JoCl  ~  r2)  RR1(a0(t),a2(t),r')dr  =  o  (3.10) 

fr2 (1  -r2)  RR1(a0(t'),a2(t),r)dr  =  0  (3.11) 

When  equation  (3.7)  is  substituted  into  (3.10)  and  (3.11)  respectively,  after  evaluating  and  simplifying  in  full  we 
obtain  the  following  systems 

_  +  JL  (_  5*0 1  +  2a0(t)  _  2 a2(t)\  8d0(t)  8d2(t)  ^  /2d0(Q  _  2a2(t)\  _  /_  _  32V021a0(t)  _ 

3  ~K£7  V  3  3  15  15  105  1  V  3  15  )  V  5  5 

16V012a2t35-  8V01a02t5-  362V01a22tl05- 278V01a0ta2tl05+16a02ta2t35+56a0ta22tl05+1475a23t3465 
+8a03tl  5+M1-  2V01 1 5-  8a0tl5-  8a2tl05=0 

(3.12) 

and 
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_  +  _±_  (_  iYo  +  2°o(t)  +  2«2(t)\  +  8d0(t)  8gz(t)  /2a0(t)  2dz(t)\  _  16^  /_  ev£  54V02a2(t) 

15  _r  RE/ V  15  15  21  /  105  315  1  V  15  21  /  V  35  35 

16V02a2t35-24V0a2t35-1107V0a2tll55-58V0a0ta2t35-334a02ta2t35-896a0ta22t385+33836a23t5005+ 
8a03t35-Ml- 2V035- 8a0tl05-  3a2t315=0. 

(3.13) 

Where  denotes  differentiation  with  respect  to  t. 

Equations  (3.12)  and  (3.13)  can  be  re-written  respectively  as 

(1848  +  9240/21)d0(t)  +  (264  -  1848/21)d2(t)  +  29568/2^ (t)  +  23600 /2af(t)  -  88704/2E01ag(t)  - 
191136/2E01a|(t)  +  29568/2a0 (t)ag(t)  +  25344/2a2(t)ag(t)  -  14678412E01a0(t)a^(t)  - 
(l848 M1  +  354816/21^  -  — )  a0(t )  -  (264MX  -  2  5  344/2 1/^  -  — )  a2(t)  =  2310G!  +  2I22Yoi  +  22176/2^  + 

v  REI  J  \  REI  J  REI 

462 MtV01  (3.14) 

and 


(3432  +  24024/21)d0(t)  +  (1144  +  17160/21)d2(t)  +  164736/2^(0  +  4872384/2af  (0  - 
494208/2 E01ao(t)  -  690768/2 V01afct)  -  6877728/2ag(t)a2(t)  -  1677312/2a|(t)a0(t)  - 

1 194336/2 E01a0(t)a2(t)  -  (3432 M1  -  1111968/2V02!  +  a0(t)  -  (429 M1  -  329472/21/q2!  +  a2(t)  = 

v  REI  J  \  REI  J 

6006G!  +  +  123552/2 Vq2!  +  2574M1701  (3.15) 

By  substituting  the  appropriate  values  of  the  parameters  Gl9  V01,  REI ,  12,  M1  and  t  into  (3.14)  and  (3.15),  after 

simplifying  we  obtain 

— 85.33333333ao(t)  -  14.77633478af  (t)  +  64.000000ag(O  +  17.0666666af(0  - 
21.33333333ag(t)a2(t)  -  85.33333333a|(t)a0(t)  +  36.57142858(t)a2(t)  -  19.14962963a0(0  - 
4.005502645a2(t)  -  2.666666667 d0(t)  +  0.5333333333d2(O  =  -3.062407407  (3.16) 

and 


— 36.57142857ao(t)  -  15.00366300ai(t)  +  27.42857142ag(t)  +  13.57575758al(t)  - 
73.14285714ag(t)a2(t)  -  54.30303030al(t)a0(t)  +  36.57142858a0(t)a2(t)  -  7.476402115a0(t)  - 
5.003597883a2(t)  -  0.5333333333 d0(t)  +  0.3809523869d2(t)  =  -0.9145767196  (3.17) 

Solving  the  system  of  non-linear  first  order  ordinary  differential  equations  (3.16)  and  (3.17)  with  initial  conditions 
a0(0)  =  0’  a2(0)  =  0  and  time  t  =  0.5  using  fourth  order  Runge-Kutta  method  we  obtain  the  values  fora0(t)  and  a2(t) 
and  when  submitted  into  (3.5)  we  obtain  the  velocity  profile  for  the  unsteady  fluid  flow  of  blood  as  shown  in  Tablet  for 
the  various  values  of  the  parameters. 

Similarly,  to  obtain  the  temperature  profile  for  the  unsteady  heat  transfer  using  the  Gerlakin’s  method,  we  assume 
a  solution  of  the  form 

0(y,  t)  =  a0(t )  +  a3(t)y  +  a4(t)y2  (3.18) 

Subjecting  (3.18)  to  the  slip  conditions  (2.15)  and  after  simplifying  we  obtain 
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6(y,t)  =  a3(t)(l-^  + a4(t)^(l-^J  (3.19) 

By  using  the  transformation  (3.3)  and  dropping  bar,  equation  (3.19)  can  be  written  as 

6(r,  t)  =  a3(t)(  1  -  r2)  +  a4(t)r2(l  -  r2)  (3.20) 


The  residual  for  the  unsteady  heat  transfer  can  be  written  as 


RR2{r,  a3(t),a4(t)) 


(3.21) 


By  substituting  (3.5)  and  (3.20)  into  (3.21),  after  simplifying  we  obtain 

/?/?2 (r,  a3 (t),  a4(t))  =  a3(t)(l  -r2)  +  a4(t)r2(l  -  r2)  -  Enl(4V^r2  -  8a0(t)V01r2  +  8a2(t)V01r2  - 
16V01a2tr2+4a02tr2-  8a0ta2tr2+16a0ta2tr4+4a22tr2-  16a22tr4+16a22tr6-  p64V01 3a2tr2-  64V013a0tr4 
-  64V01 3a2tr4-  128V013a2tr6+  96V012a22tr4 + 96V012a02tr4-  384V012a22tr6+384V012a22tr8-  512V01a23 
trl 0-  384V01a23tr6-  64V01a03tr4+768V01a23tr8+64V01a23tr4-  64a03ta2tr4+128a03ta2tr6+96a02ta22t 
r4-  384a02ta22tr6 +384a02ta22tr8+512a0ta23trl 0+384a0ta23tr6-  768a0ta23tr8-  64a0ta23tr4+16a23tr4 
- 1 28a23r6 + 1 6a04r4-  512a24rl0+384a24r8+16V014r4+256a24rl2+192V01a02a2r4-  384V01a02a2r6- 1 92 
V 01  a0a22r4 + 768V 01  aO a2 2 r 6-  768V01a0a22r8- 1 92V012a0a2r4-  (pi-  4V01ra0t+4V01ra2t-  8r4V01a2t+4r2 
aOtaOt-  4r2a0ta2t+8r4a0ta2t-  4r2a2ta0t+4r2a2ta2t-  8r4a2ta2t+8r4a2ta0t-  8r4a2ta2t+l  6r6a2ta2t-Al  4a 
3t-  4a4t+l  6a4tr2  (3 .22) 


By  replacing  ft/?!(a0(t),  a2(t),r)  with  RR2(a3(t),  a4(t),  r)  in  (3.10)  and  (3.11),  after  evaluating  and  simplifying 
in  full  we  obtained  the  following  system  of  equations 


8a3(t) 

15 


+  2fiM  +  4£ni 
105  nl 


(  2V$i 
V  15 


+  4Vb1fl0(t)  60  V01a2(t) 


15 


105 


-2ap(t)  4flo(t)a2(t) 
15  105 


64V031a2(t) 

315 


8V013a0t35- 228V012a22t3465-  12V012a02t35+504V01a23t45045+8V01a23t35-  8a2ta03t315-  56a0ta23t5 
005-228a02ta22t3465-2a04t35-3090a24t675675+8p0a2ta02tl05+152V01a0ta22tll55+24V0124a0ta2t35 


+4  pi  V01a0t4-  51  V01a2tl  40-  6a2ta2tl89+2a0ta2tl05-  68a2ta0tl05+4A-  2a3t3+2a4tl5=0 

(3.23) 


and 


803(0 

105 


+  2£l “  +  4Fnl 

315  nl 


(  2Vqi 
V  35 


|  4Vb1fl0(t)  4Vb1a2(Q 


35 


315 


2 apCQ  4aQ(t)a2(t) 

315  315 


38q|(Q\  16  ,  (  64^02(0 

3465  )  ^  V  V  693 


8J7013a0t63- 372V012a22t9009-  12V012a02t63+l 992V01a23tl35135+8V01a23t63-  8a2ta03t231- 1 992a0ta 


231135135- 1 24a02ta22t3 003- 2a04t63-  7074a24t297295+22p0a2ta02t6935+744V01a0ta22t9009+24V012 


4a0ta2t63+4pl  V01a0tl2- 222V01a2tl512-  86a2ta2tll55+2a0ta2t315-  24a2ta0t315+4A-  2a3tl5- 1 0a4tl  0 
5=0  (3.24) 


Equations  (3.23)  and  (3.24)  can  be  re-written  respectively  as 

6756750Voido(t)  -  984555 4>V01d2(t)  +  5148O0Void2(t)  -  17503200!  a2(t)a0(t)  - 
858OO01a2(t)a2(t)  +  61776O0a£(t)  -  2O8783575O0a|(t)  +  2471O4O0l7ola^(t)  +  120960 0Vrola|(t)  - 
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120960 0%(t)af(t)  -  274560 (pa2{t)al{t)  +  1422720 $K01a0(t)a|(t)  +  823680 (pV01a2{t)a20{t)  - 
71136O0aj;(t)a|(t)  -  (360360 Enl  +  37O656O0Vo21)aj;(t)  -  (85800£nl  +  711360 0ko21)a|(t)  +  ( 720720EnlV01  + 
2471040 0i7ko31)%(t)  +  (1544400 EnlV01  -  10810800 ^K03i)a2(t)  +  (1544400Fnl  +  741312O0To21)%(t)a2(t)  + 
(360360  -  18018004)%  (t)  +  (51480  +  3603604)a4(t)  =  360360£nlF021  (3.25) 


and 


765765 (f>V01a0(t)  —  1349205 04%(t)  —  5834401%(t)%(t)  +  3734O1601%(t)%(t)  — 

684216 01%(t)a2(t)  -  1166880 ^a^(t)  -  1131840a|(t)  +  4667520 (pvV01a%(t)  +  5418240Folaf  (t)  - 
5418240%(t)af  (t)  —  127296O0%(t)ao(t)  +  116688O0Fol%(t)ao(t)  +  3O3552O0Fol%(t)a|(t)  — 
151776O0a^(t)a|(t)  -  (525096Fnl  +  700128001^)%;  (t)  -  (100776Enl  +  151776O0Fo2L)af  (t)  + 
(1050192£'nlV01  +  466752O0Vo31)ao(O  +  (116688EnlV01  -  339456O0To31)%(t)  - 
(116688Fnl  +  14OO256O07o2i)ao(Oa2(O  +  (175032  -  12252244)a3(0  +  (58344  -  8751604)a4(0  = 
525096£'nlV021  (3.26) 

Substituting  the  appropriate  values  of  the  parameters  0,  01;  V01,  Enl  and  Ax  into  (3.25)  and  (3.26)  respectively, 
we  obtain 

0. 1066666667% (t)  -  0.01523809524d2(O  +  0.53333333%(t)  +  0.0761904769d4(0  + 

0.076190476%  (t)%(t)  +  0.076190476%(t)%(t)  +  0.53333333%(t)%(t)  -  0.126984127d2(0a2(0  - 
54.857142857(4(1)  -  4.3902763903a|(O  +  43.88571429<4(t)  +  2.1148251748a|(t)  + 

24.38097381ao(t)%(t)  -  10.741258741%(0a|(0  -  13.96571429(4(1)  +  14.62857143ag(0ag(0  + 

2 5.26753247% (t)a|(t)  -  2.697142857%(0%(0  +  4.9333333%(0  -  0.98666666%(0  =  0.1197714286  (3.27) 

and 


0. 04571428571%  (t)  +  0.005079365099%(t)  +  0.076190476%(t)  +  0.0253968254d4(O  - 
0.0253968254%(t)%(t)  -  0.0253968254%(t)%(t)  -  0.2285714286%(t)%(t)  -  0.0438672439%(t)%(t)  - 
30.476190476(4(1)  -  2.9561027208(4(1)  +  24.38095238(4(0  +  2.830236430(4(1)  - 
33.246753247ag(0%(0  -  14.151182151%(Oaf  (0  -  7.657142857(4(0  -  1.651415252(4(1)  + 
19.94805195ag(0%(0  +  15.85614386%(t)af  (t)  -  39.64035964(4  (0<4(t)  +  1.112380952%(O  + 

0.2812 12 121 2% (t)  -  4.065800866% (t)% (t)  +  0.98666666%(0  +  0.7047619048%(0 

=  0.06247619048  (3.28) 

Using  fourth  order  Runge-kutta  method,  we  solve  the  system  of  non-linear  ordinary  differential  equations  (3.16), 
(3.17),  (3.27)  and  (3.28)  with  initial  conditions  %(0)  =  0, %(0)  =  0, %(0)  =  0,%(0)=0  and  time  1  =  0.  When 
substituting  the  values  of  %(t)  and  %(t)  into  (3.20)  after  simplifying  we  obtain  the  temperature  profile  for  the  unsteady 
heat  transfer  as  shown  in  Table  2  for  various  values  of  the  parameters. 

Volume  Flow  Rate 

The  volume  flow  rate  denoted  by  Q  is  given  by 

Q  =  2n  rw(r)dr  (3.29) 


Impact  Factor  (JCC):  4.9784 
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Putting  (4.7)  into  (3.44)  and  evaluate  to  obtain 

Q  =  12  [3V0(R(z))4  +  a0(t )  (6(RO))2  -  3(/?(z))4)  +  a2(t)  (3(/?0))4  -  2 (/?(z))6)]  (3.30) 

Shear  Stress 


The  shear  stress  denoted  by  Tsis  given  as 


r=R(z ) 


+ 


:f  I 

/  I  r=R(z ) 


(3.31) 


Simplified  (3.46)  to  obtain 

rs  =  2/i/?(Z)  (V0  -  a0(t)  +  a2(t)  -  2R((Z))2a2(t)}  +  16R(Z)/33(K0  -  a0(t)  +  a2(t)  -  2 (R(Z))2a2)  (3.32) 

Resistance  to  Flow 


The  resistance  to  flow  can  be  denoted  as  x\i  and  is  given  by 

_dP_ 

\l/  =  - — - 

1 2  [3  V0  {R  (z))4 + a0  (t)  (6  (i?  (z))2  -  3  (ft  (z))4) + a2  (t)  (3  (i?  (z))4  -  2  (7?  (z))  6)] 


(3.33) 


Table  1:  Values  of  the  Parameters  used  in  the  Numerical  Results  and  the  Corresponding  Velocity 


Profile  for  the  Unsteady  Blood  Flow  Model 


Figures 

Gi 

Voi 

REI 

n 

il. 

t 

w(r,  t) 

1.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.1437r2  +  0.3437  +0.0335r2(l-r2) 

2 

2.0 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.1858r2  +  0.3858  +0.0296r2(l-r2) 

2.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.2290r2  +  0.4290  +0.0262rz(l-r2) 

1.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.2261r2  +  0.4261  +0.0260r2(l-r2) 

3 

1.5 

0.25 

0.9 

10 

1 

0.65 

0.5 

-0.2174r2  +  0.4174  +0.0249r2(l-r2) 

1.5 

0.25 

0.9 

10 

1 

0.95 

0.5 

-0.2090r2  +  0.4090  +0.0237r2(l-r2) 

1.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.0572r2  +  0.2572  +0.0009r2(l-r2) 

4 

1.5 

0.35 

0.9 

10 

1 

0.35 

0.5 

-0.1538r2  +  0.3538  +0.0019r2(l-r2) 

1.5 

0.45 

0.9 

10 

1 

0.35 

0.5 

-0.25 10r2  +  0.4510  +0.0030r2(l-r2) 

1.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.1246r2  +  0.3246  +0.0168r2(l-r2) 

5 

1.5 

0.25 

0.9 

10 

5 

0.35 

0.5 

-0.20 16r2  +  0.4016  +0.0073r2(l-r2) 

1.5 

0.25 

0.9 

10 

9 

0.35 

0.5 

-0.2542r2  +  0.4542  +0.0001r2(l-r2) 

1.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.3246r2  +  0.5246  +0.0499r2(l-r2) 

6 

1.5 

0.25 

0.9 

20 

1 

0.35 

0.5 

-0.3108r2  +  0.5108  +0.0538r2(l-r2) 

1.5 

0.25 

0.9 

30 

1 

0.35 

0.5 

-0.2848r2  +  0.4848  +0.0554r2(l-r2) 

1.5 

0.25 

0.3 

10 

1 

0.35 

0.5 

-0.0053r2  +  0.2053  +0.000Ir2(l-r2) 

7 

1.5 

0.25 

0.6 

10 

1 

0.35 

0.5 

-0.0094r2  +  0.2094  +0.0001r2(l-r2) 

1.5 

0.25 

0.9 

10 

1 

0.35 

0.5 

-0.009  lr2  +  0.2091  +0.0005r2(l-r2) 

Table  2:  Values  of  the  Parameters  used  in  the  Numerical  Results  and  the 
Corresponding  Temperature  Profile  for  the  Unsteady  Heat  Transfer  Model 


Figures 

<!> 

Enl 

A, 

^01 

9(r,t) 

1.25 

1.5 

1.35 

0.25 

2.0 

0.007  -  0.007r2+  0.001r2(l-r2) 

8 

1.50 

1.5 

1.35 

0.25 

2.0 

0.008  -  0.008r2  +  0.013r2(l-r2) 

1.75 

1.5 

1.35 

0.25 

2.0 

0.009  -  0.009r2  +  0.015r2(l-r2) 

1.25 

1.5 

1.35 

0.25 

2.0 

0.006  -  0.006r2  +  0.008r2(l-r2) 

9 

1.25 

1.8 

1.35 

0.25 

2.0 

0.007  -  0.007r2+  0.009r2(l-r2) 

1.25 

2.1 

1.35 

0.25 

2.0 

0.008  -  0.008r2+  0.010r2(l-r2) 

1.25 

1.5 

1.35 

0.25 

2.0 

0.007  -  0.007r2+  0.008r2(l-r2) 

10 

1.25 

1.5 

1.65 

0.25 

2.0 

0.006-  0.006r2+  0.006r2(l-r2) 

1.25 

1.5 

1.95 

0.25 

2.0 

0.005  -  0.005r2+  0.005r2(l-r2) 
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Table  2  Contd., 

1.25 

1.5 

1.35 

0.25 

1.0 

0.006  -  0.006r2  +  0.008r(  1  -r2) 

11 

1.25 

1.5 

1.35 

0.25 

2.0 

0.010-  O.OlOr2  +  0.001r2(l-r2) 

1.25 

1.5 

1.35 

0.25 

3.0 

0.013  -  0.013r2+  0.012r2(l-r2) 

1.25 

1.5 

1.35 

0.25 

2.0 

0.005  -  0.005r2+  0.006r2(l-  r2)  0.006  - 

12 

1.25 

1.5 

1.35 

0.35 

2.0 

0.006r2  +  0.007r2(l-r2) 

1.25 

1.5 

1.35 

0.45 

2.0 

0.007  -  0.007r2+  0.008r2(l-r2) 

Figure  2:  Variation  of  Velocity  Profile  of  the  Unsteady  Blood  Flow  Model 
with  Increasing  Values  of  the  Pressure  Gradient  in  the  Radial  Direction. 


Figure  3:  Variation  of  Velocity  Profile  of  the  Unsteady  Blood  Flow  Model  with 
Increasing  Values  of  the  Magnetic  Field  Parameter  in  the  Radial  Direction. 


Impact  Factor  (JCC):  4.9784 


NAAS  Rating  3.45 
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radial  distance^) 


- ^or'0-25 - vor°3S  —  For0-45 


Figure  4:  Variation  of  Velocity  Profile  of  the  Unsteady  Blood  Flow  Model 
with  Increasing  Values  of  the  Slip  Velocity  in  the  Radial  Direction. 


Figure  5:  Variation  of  Velocity  Profile  of  the  Unsteady  Blood  Flow  Model 
with  Increasing  Values  of  the  Shear  Thickening  in  the  Radial  Direction. 
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radial  distance^} 


- 0=10 - 0=20 - 0=30 


Figure  6:  Variation  of  Velocity  Profile  of  the  Unsteady  Blood  Flow  Model 
with  Increasing  Values  of  the  Shear  Thinning  in  the  Radial  Direction. 


radial  distance^) 

-  RE  1=0.3 - RE  1  =0 .6 - RE  1=0 .9 


Figure  7:  Variation  of  Velocity  Profile  of  the  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Reynolds  Number  in  the  Radial 

Direction. 


Impact  Factor  (JCC):  4.9784 
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radial  distance^) 


0=1.25  - 0=1 .50 -  0=1.75 


Figure  8:  Variation  of  Temperature  Profile  of  the  Unsteady  Heat  Transfer 
with  Increasing  Values  of  the  Shear  Thinning  in  the  Radial  Direction. 


Figure  9:  Variation  of  Temperature  Profile  of  the  Unsteady  Heat  Transfer 
with  Increasing  Values  of  the  Eckert  Number  in  the  Radial  Direction. 
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Figure  10:  Variation  of  Temperature  Profile  of  the  Unsteady  Heat  Transfer 
with  Increasing  Values  of  the  Third  Grade  Parameter  in  the  Radial  Direction. 


Figure  11:  Variation  of  Temperature  Profile  of  the  Unsteady  Heat  Transfer 
without  Hematocrit  with  Increasing  Values  of  the  Shear  Thickening  in  the  Radial 

Direction. 


Impact  Factor  (JCC):  4.9784 
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Figure  12:  Variation  of  Temperature  Profile  of  the  Unsteady  Heat  Transfer  with 
Increasing  Values  of  the  Slip  Velocity  in  the  Radial  Direction. 


Figure  13:  Variation  of  Temperature  Profile  of  the  Unsteady  Heat  Transfer  with 
Increasing  Values  of  the  Magnetic  Field  Parameter  in  the  Radial  Direction. 
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Figure  14:  Variation  of  Velocity  Profile  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  Time  along  the  Radial  Direction. 


Figure  15:  Variation  of  Temperature  Profile  of  Unsteady  Heat  Transfer 
Model  Flow  Model  with  Increasing  Values  of  Time  along  the  Radial  Direction. 


Impact  Factor  (JCC):  4.9784 
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- =0.001 - =0.00  2  - F01=0.003 


Figure  16:  Variation  of  Volume  Flow  Rate  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Slip  Velocity  in  the  Entire 
Stenotic  Region  along  the  Axial  Direction. 


F01  =  0.001 - Vm  =  0.002  - Vm  =  0.003 


Figure  17:  Variation  of  Wall  Shear  Stress  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Slip  Velocity  in  the  Entire 
Stenotic  Region  along  the  Axial  Direction. 
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F01=0.001 - F01-0.002  - ^=0.003 


Figure  18:  Variation  of  Wall  Shear  Stress  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Slip  Velocity  in  the  Entire 
Stenotic  Region  along  the  Axial  Direction. 


3^=0.005  - 3^=0.010 - A^  =0.0 15 


Figure  19:  Variation  of  Volume  Flow  Rate  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Magnetic  Field  Parameter  in  the 
Entire  Stenotic  Region  along  the  Axial  Direction. 


Impact  Factor  (JCC):  4.9784 


NAAS  Rating  3.45 
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Mj  =0.005  - =0.010 - A^  =0.0 15 


Figure  20:  Variation  of  Wall  Shear  Stress  of  Unsteady  Blood  Flow  Model 
with  Increasing  Values  of  the  Magnetic  Field  Parameter  in  the  entire 
Stenotic  Region  along  the  Axial  Direction. 


A/|=0 .00  5  - A/|=0 .010 - M{= 0.015 


Figure  21:  Variation  of  Resistance  to  Unsteady  Blood  Flow  Model 
with  Increasing  Values  of  the  Magnetic  Field  Parameter  in  the 
Entire  Stenotic  Region  along  the  Axial  Direction. 


29 


www.iaset.us 


editor@iaset.us 


30 


Jimoh  A,  Okedayo  G.T  &  Ikpakyegh  L.N 


Q=2.0 - 0=2.05 - 0=2.10 


Figure  22:  Variation  of  Volumetric  Flow  Rate  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Shear  Thinning  in  the  Entire 
Stenotic  Region  along  the  Axial  Direction. 


0=2.0 - 0=2.05 - 0=2.10 


Figure  23:  Variation  of  Wall  Shear  Stress  of  Unsteady  Blood 
Flow  Model  with  Increasing  Values  of  the  Shear  Thinning  in 
the  Entire  Stenotic  Region  along  the  Axial  Direction. 


Impact  Factor  (JCC):  4.9784 


NAAS  Rating  3.45 
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- 0=2.0 - 0=2.05 - 0=2.10 


Figure  24:  Variation  of  Resistance  to  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Shear  Thinning  in  the 
Entire  Stenotic  Region  along  the  Axial  Direction. 


- Q=l.Q - Qj=1.02 - Qj=1.04 


www.iaset.us 


Figure  25:  Variation  of  Volumetric  Flow  Rate  of  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Shear  Thickening  in  the  Entire 
Stenotic  Region  along  the  Axial  direction. 
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Figure  26:  Variation  of  Wall  Shear  Stress  of  Unsteady  Blood 
Flow  Model  with  Increasing  Values  of  the  Shear  Thickening  in 
the  entire  Stenotic  region  along  the  Axial  Direction. 
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Figure  27:  Variation  of  Resistance  to  Unsteady  Blood  Flow 
Model  with  Increasing  Values  of  the  Shear  Thickening  in  the 
Entire  Stenotic  Region  along  the  Axial  Direction. 


DISCUSSIONS  OF  RESULTS 

In  this  work,  we  have  studied  the  effects  of  velocity,  magnetic  field  and  other  parameters  on  the  flow  and  heat  transfer 
characteristics  by  considering  the  blood  as  a  third  grade  fluid,  and  the  results  based  on  the  mathematical  analysis  indicate 
that  an  increase  in  slip  velocity  significantly  leads  to  an  increase  in  flow  velocity,  flow  rate  and  shear  stress  as  shown  in 
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Figures  4,  16,  17  respectively.  In  the  same  vein,  heat  transfer  and  resistance  to  flow  decreases  as  the  slip  velocity  increases, 
these  are  depicting  in  figures  12  and  18  respectively.  Increase  in  magnetic  field  parameter  leads  to  decease  in  flow 
velocity,  temperature  profile,  flow  rate,  shear  stress  but  increase  the  resistances  to  fluid  flow  as  indicated  in  Figures  3,  13, 
19,  20  and  21,  respectively.  Pressure  gradient  and  Reynolds  number  increase  with  the  flow  velocity  as  shown  in  Figures  2 
and  7,  respectively.  Eckert  number  increase  with  the  temperature  profile  as  indicated  in  figure  9  while  increase  in  third 
grade  parameter  leads  to  decrease  in  temperature  profile  as  shown  in  figure  10. 

Other  important  properties  of  blood  are  shear  thinning  and  shear  thickening.  Increase  in  shear  thinning  leads  to 
increase  in  the  flow  velocity,  temperature  profile,  flow  rate  and  shear  stress  but  reduces  the  flow  resistance  as  shown  in 
figures  6,  8,  22,  23,  and  24,  respectively.  Flow  velocity,  flow  rate  and  shear  stress  decrease  as  shear  thickening  increases 
and  these  are  shown  in  figures  5,  25,  and  26,  respectively  while,  temperature  profile  and  resistance  to  flow  increase  with 
shear  thickening  as  shown  in  figures  11  and  27,  respectively.  Finally,  velocity  and  temperature  profiles  increase  with  time 
as  indicated  in  figures  14  and  15,  respectively. 

CONCLUSIONS 

The  computational  analysis  of  unsteady  flow  of  blood  and  heat  transfer  through  an  artery  with  stenosis  and  slip  conditions 
have  been  studied  through  this  paper.  The  mathematical  problems  are  solved  analytically  and  the  significant  findings  are 
summarized  below: 

•  Velocity  profile,  volumetric  flow  rate,  and  shear  stress  increase  while  the  temperatures  profile  and  resistance  to 
flow  decreases  with  increasing  values  of  the  slip  velocity. 

•  Velocity  profile,  temperature  profile,  volumetric  flow  rate,  and  shear  stress  decrease  while  resistance  to  flow 
increase  as  magnetic  field  parameter  increases. 

•  Velocity  profile,  temperature  profile,  flow  rate,  and  shear  stress  increase  while  resistance  to  flow  decreases  with 
increasing  values  of  the  shear  thinning. 

•  Velocity  profile,  flow  rate,  and  shear  stress  decrease  while  temperature  profile  and  resistance  to  flow  increase  as 
the  values  of  the  shear  thickening  increases. 

•  Finally,  increase  in  time  positively  influence  both  the  flow  and  heat  transfer  rate. 

The  research  analysis  incorporating  externally  applied  magnetic  field  is  useful  for  the  reduction  of  blood  flow 
during  surgery  and  magnetic  resonance  imaging  (MRI).  Also,  the  incorporating  slip  velocity  in  the  constricted  artery  can 
help  to  reduce  the  blood  viscosity  which  in  effect  influences  the  flow  of  blood.  It  is  clear  that  the  magnetic  field  and  slip 
velocity  are  the  strong  parameters  influencing  the  flow.  Atherosclerosis  leads  to  hypertension  because  the  deposition  of 
cholesterol  or  micro  molecules  at  the  arterial  wall  will  not  only  reduce  the  regular  blood  supply  but  also  cause  the 
hardening  or  thickening  of  the  blood  vessel  in  a  way  that  the  blood  vessel  cannot  contract  easily.  In  view  of  these,  the 
present  study  can  also  be  useful  to  control  the  blood  flow  in  disease  state. 
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NOMENCLATURES 

w-Fluid  velocity 
w-Dimensionless  fluid  velocity 
t-Time  component 
F-Dimensionless  time  component 
r-Radial  distance 
y-Dimensionless  radial  distance 
z- Axial  distance 
ws-Slip  velocity 

VQ1 -Dimensionless  Slip  velocity  for  the  unsteady  blood  flow 

T-Temperature  profile 

7W-Pipe  temperature 

0  -Dimensionless  temperature  profile 

Tm -Fluid  temperature 

R0  -Radius  of  the  normal  artery 

/?o -Magnetic  Field  Strength 

R(z)-Radius  of  the  artery  in  a  stenotic  region 

cr-Electrical  Conductivity 

^-Resistance  to  flow 

K-Thermal  conductivity 

Q- Volume  trie  flow  rate 

rs-Wall  Shear  Stress 
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^-Maximum  height  of  the  stenosis 
L-Length  of  the  stenosis 
W-Fluid  velocity 

Gr -Pressure  gradient  for  the  unsteady  flow 

Mr -Magnetic  field  parameter  for  the  unsteady  flow  Q-Shear  thinning  for  the  unsteady  flow  -Shear  thickening  for  the 
unsteady  flow 

Enl  -Eckert  number  for  the  unsteady  heat  transfer 
c[)-Shear  thinning  for  the  unsteady  heat  transfer 
§  ■ -Shear  thickening  for  the  unsteady  heat  transfer 
Ax -Third  grade  parameter  for  the  unsteady  heat  transfer 
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